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Purpose: Retinoblastoma (RB) sets the paradigm for hereditary cancer syndromes, for which medical care can change 
depending on the results of genetic testing. In this study, we screened constitutional mutations in the RBI gene via a 
method combining DNA sequencing and multiplex ligation-dependent probe amplification (MLPA), and performed a 
preliminary exploration of genotype-phenotype correlations. 

Methods: The peripheral blood of 85 retinoblastoma probands, including 39 bilateral and 46 unilateral, was collected, 
and genomic DNA was extracted. DNA sequencing was conducted first. MLPA analysis was applied for patients with 
bilateral RB with negative sequencing results and unilateral probands whose age at diagnosis was less than 1 year old. 
Results: Thirty-four distinct mutations were identified in 40 (47.1%) of the 85 probands (36 bilateral and four unilateral), 
of which 20% (8/40) was identified by MLPA. The total detection rate in bilateral cases was 92.3% (36/39). Of the total 
mutations identified, 77.5% (31/40) probands with a mean age of 10.7 months at diagnosis had null mutations, and 22.5% 
(9/40) with a mean age of 13.5 months at diagnosis had in-frame mutations. Of the 31 probands with null mutations, 
bilateral RB accounted for 96.8% (30/31). Of the nine probands with in-frame mutations, 66.7% had bilateral RB. There 
were seven new mutations of RBI identified in this report, including six null mutations and one missense mutation. Clini- 
cal staging of the tumor did not show obvious differences between patients with null mutations and in-frame mutations. 
Conclusions: Our results confirm that the type of mutation is related to age of onset and the laterality, but not staging 
of the retinoblastoma tumor. MLPA is a reliable method for detecting gross deletion or duplication of the RBI gene. The 
combination of sequencing and MLPA improves the clinical diagnosis of RB. 



Retinoblastoma (RB; OMIM 180200) is the most 
frequent primary intraocular malignant tumor in children, 
probably arising from cone precursor cells [1]. RB mainly 
affects children under 6 years old, with an incidence rate of 1 
case per 15,000 to 20,000 live births [2,3]. According to the 
"two-hit" hypothesis, RB includes hereditary and nonheredi- 
tary forms, resulting from the mutation of both alleles of the 
RBI gene (Gene ID: 5925) [4,5]. Approximately 40% of the 
cases (including all bilateral and 15% of unilateral RB) are 
heritable, carrying a germ-line RBI mutation transmitted as 
an autosomal dominant trait with 90% penetrance [6,7]. The 
other 60% of cases are non-heritable RB (85% unilateral RB), 
caused by the inactivation of both alleles in the developing 
retina [8]. 
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The RBI gene was the first tumor-suppressor gene 
found, located in 13ql4.2, and the whole DNA length is 183 
kb. Mutations in the RBI gene are highly heterogeneous 
and scattered in the promoter and the 27 coding exons. To 
date, more than 1,600 distinct mutations, ranging from small 
mutations to large deletions, have been registered in the RBI 
Gene Mutation Database. In developed countries, RBI gene 
testing has already been applied as a routine examination 
in RB probands [9,10], since the detection of RBI mutations 
provides evidence for genetic counseling and clinical manage- 
ment. However, molecular diagnosis is still in its early stages 
in China. A search of the PubMed and Chinese databases, 
including CNKI (China National Knowledge Infrastructure) 
and WanFang data, revealed several studies regarding RBI 
alterations in Chinese patients with RB [11-21]. However, 
these studies mainly focus on pedigree or have a small sample 
size, with no definitive analysis of mutational characteris- 
tics and the correlation between genotype and phenotype. 
Therefore, our cooperative group, which is based on encour- 
aging treatment outcomes of patients with RB reported in a 
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previous article [6], developed the molecular detection project 
to further enhance the level of clinical management of retino- 
blastoma. We have already reported several RBI alterations 
by using DNA sequencing with a detection rate of 78.6% in 
bilateral RB [22], which had not been capable of screening 
all RBI variations. Here, we demonstrate that combining 
DNA sequencing and multiplex ligation-dependent probe 
amplification (MLPA) for detecting the mutation spectrum 
of RBI allows for the preliminary exploration of genotype- 
phenotype correlations. 

METHODS 

Patients: This study recruited a total of 85 unrelated RB 
probands, including 37 (43.5%) boys and 48 (56.5%) girls, 
who were diagnosed between January 2012 and October 2013 
in a Chinese cooperative group consisting of the Children's 
Hospital of Fudan University and the Eye & ENT Hospital 
of Fudan University. The genetic testing results, using 
Sanger sequencing, of 35 of these patients were reported in 
a previous paper published in Chinese [22], but new find- 
ings emerged using the complementary MLPA method. The 
combined results are discussed in the present paper. The 
patients included 39 cases of bilateral RB and 46 cases of 
unilateral RB. Clinical features, including laterality, age at 
diagnosis, and International Classification of Intraocular 
Retinoblastoma (ICRB) staging at diagnosis (Appendix 1), 
were recorded. After the study protocol was approved by the 
Institutional Review Board of Fudan University, informed 
consent for genetic testing was obtained from the guardians 
of children affected with RB. The study protocol was in 
accordance with the provisions of the Declaration of Helsinki 
and the ARVO (the Association for Research in Vision and 
Ophthalmology) statement on human subjects. Peripheral 
blood samples (3 ml) were collected from the probands in 
EDTA (EDTA) anticoagulant tubes. Blood samples were 
stored at -20 °C until DNA was extracted. 

RBI gene sequencing: Genomic DNA was isolated from 
peripheral blood leukocytes using the QIAamp DNA Blood 
Mini Kit (Qiagen, Dusseldorf, Germany), according to the 
manufacturer's instructions. Extracted DNA was measured 
with spectrophotometry (NanoDrop 2000; Thermo Scientific, 
Wilmington, DE) to ensure a concentration of 40 ng/(a.l. 

Primers were designed according to a previous study by 
Abouzeid et al., with some modifications, using the Oligo 
software (National Biosciences, Plymouth, MN) [23]. Twenty- 
four pairs of primers were generated, which covered all 
coding exons and at least 50 bp flanking intronic sequences 
of each exon. PCR reactions were performed in a thermal 
cycler (ABI 2700), in a total volume of 10 (j.1 containing 40 ng 
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of genomic DNA, 2 picomoles of each primer, 0.2 \il of Taq 
DNA polymerase, and 1 (il of 10 x buffer. Dimethyl sulfoxide 
(DMSO, 10%) was added to the exon 1 amplification system. 
Reactions were performed for the first 14 cycles of denatur- 
ation at 95 °C for 30 s, annealing at a specific temperature 
gradient (65 °C, -0.5 °C per cycle) for 30 s, extension at 72 °C 
for 40s; the second 20 cycles of denaturation at 95 °C for 30 
s, annealing at 58 °C for 30 s, extension at 72 °C for 40 s, and 
a final extension step at 72 °C for 5 min. The PCR products 
were qualified with agarose gel electrophoresis. 

After the unincorporated deoxynucleoside triphos- 
phates and primers were removed using Exonuclease I 
and shrimp alkaline phosphatase, the PCR products were 
directly sequenced using the ABI BigDye Terminator v 3.1 
Sequencing Standard Kit (Carlsbad, CA) and run on an ABI 
3500 Genetic Analyzer. The sequence data were analyzed 
by comparison to the consensus sequence of the RBI gene 
(NM_000321), using Mutation Surveyor (V 4.0) software 
(SoftGenetics; State College, PA). Additional information was 
obtained from the literature, the RBI Gene Database (LOVD, 
2010), the Human Gene Mutation Database (HGMD, 2007), 
and the Thousand Genomes Project (Thousand Genome 
Consortium, 2010). The functional consequences of the 
mutations were predicted by bioinformatics tools, SIFT, and 
PolyPhen-2. 

Multiplex ligation-dependent probe amplification analysis: 
MLPA analysis was applied in bilateral probands who did 
not show germ-line mutations in DNA sequencing and in 
unilateral probands whose guardians demanded additional 
genetic examinations, considering the earlier age at diagnosis, 
which was less than 1 year old, although the DNA sequencing 
results had been negative. Gross deletion and duplication of 
RBI were examined using methylation-specific (MS)-MLPA 
analysis, the commercially available SALSA MLPA P047 
RBI probemix (MRC-Holland, Amsterdam, The Nether- 
lands). This P047-C1 RBI probemix contains probes for 26 
of the 27 RBI exons and 13 reference probes. This method 
was used to detect loss of heterozygosity (LOH) status of 
the RBI gene located on 13ql4. The procedure was based 
on the manufacturer's protocol. Each MS-MLPA reaction 
was subjected to further analysis on an ABI 3730 automated 
sequencer (Applied Biosystems). The data were analyzed 
with GeneMarker vl.91 (SoftGenetics) according to the 
manufacturer's instructions. The relative value of probes 
compared to the reference control was calculated. Threshold 
ratios for deletion and duplication were set at 0.75 and 1.35, 
respectively. The peak ratios for all target loci in a 1:1 ratio 
indicated normal copy numbers. The reduced or amplified 
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signal of probes with <0.75 or >1.35 indicated loss or gain, 
respectively, of the copy number in these probe loci. 

Statistical analysis: Data were analyzed using SPSS software 
(version 16.0, Chicago, IL). Statistical analysis was performed 
with one-way ANOVA, the chi-square test, and Fisher's exact 
test. p<0.05 was considered to indicate statistical significance. 

RESULTS 

Patients included 39 (45.9%) bilateral and 46 (54.1%) unilat- 
eral cases, none of whom had a family history of RB. The 
age at diagnosis ranged from 1 to 72 months, with a median 
age of 11.1 months for patients with bilateral RB and 27.1 
months for patients with unilateral RB. By combining the 
two different approaches, we discovered 34 distinct mutations 
(47.1%) in 40 probands (36 bilateral and four unilateral), of 
which 20% (8/40) was identified by MLPA. The ability to 
identify a germ-line mutation from blood samples of patients 
with bilateral RB was 92.3% (36/39). Of the total mutations 
identified (Table 1) in 40 patients, eight (20%) had gross dele- 
tions involving multiple exons or the whole gene, 17 (42.5%) 
had nonsense mutations, four (10%) had small deletions or 
insertions, two (5%) had splice mutations, and nine (22.5%) of 
patients had missense mutations, all of which were heterozy- 
gous mutations, except one homozygous missense mutation, 
c.2455C>G. The details of each group are as follows. 

Null mutation: A null mutation is a pathogenic mutational 
type, which refers to those with mRNA harboring a prema- 
ture translational-termination codon, including whole allele 
deletions, deletions involving one or more exons, and splice 
mutations, small frame-shifting intraexonic insertions and 
deletions, and nonsense mutations [24,25]. Our group identi- 
fied 25 distinct null mutations in 31 (36.5%) of the 85 patients 
(Table 1). Of these mutations, the C.17350T mutation was a 
recurrent nonsense mutation identified in three patients. Six 
null mutations had not been reported before. The 31 patients 
with RB with an average age of 10.7 months at diagnosis 
included one patient (3.2%) with unilateral RB and 30 patients 
(96.8%) with bilateral RB (Table 2). Of the 61 eyes studied in 
this group, 20 (32.8%) presented with group E ICRB followed 
by groups C, D, B, and A ICRB, and the other four (6.6%) 
presented with extraocular extension at diagnosis (Table 3). 

In-frame mutations: In-frame mutations stand for those not 
truncating based on prediction, such as indels (small inser- 
tions and deletions) in frame, and missense mutations [24]. 
We identified nine missense mutations in nine probands, 
of which one was a new mutation. Predicting the effects of 
these non-synonymous variants on protein function indicated 
'damaging' designations for six patients with bilateral RB 
and 'tolerated' for three patients with unilateral RB (Table 
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1). The nine probands with an average age of 13.5 months 
at diagnosis involved 66.7% (6/9) bilateral and 33.3% (3/9) 
unilateral patients (Table 2). Of the 15 eyes, eight (53.3%) 
presented with group E ICRB, followed by groups D, A/B, 
and C ICRB, and one (6.7%) presented with extraocular 
extension (Table 3). 

Negative genetic testing results: This group represented 
patients with negative genetic testing results. There were 
45 cases, including three (6.67%) bilateral and 42 (93.3%) 
unilateral. The mean age at diagnosis was 27.3 months. Of the 
48 eyes, 23 (47.9%) presented with group E ICRB, followed 
by groups D/C, B, and A ICRB. Sixteen of the 48 (33.3%) 
presented with extraocular extension (Table 3). 

DISCUSSION 

Genotype- phenotype correlations and implications for 
genetic counseling: In our study, all mutations were classi- 
fied based on the prediction whether the transcripts harbored 
a premature translational-termination codon. Transcripts 
from alleles of the RBI gene with null mutations have been 
shown to be subject to degradation by nonsense-mediated 
decay (NMD) [25-27], resulting in no detectable protein. 
Therefore, a null mutation is a comparatively more severe 
mutation than an in-frame mutation. Of the 31 probands with 
null mutations, 96.8% had bilateral RB, while only 66.7% 
of the probands with in-frame mutations developed bilateral 
RB (Table 2). Early onset retinoblastoma was seen more with 
null mutations (average age at diagnosis was 10.7 months) 
than in-frame mutations (average age at diagnosis was 13.5 
months). These results are consistent with reports from 
Canada and India [24,28], and confirm that null mutations 
tend to the severe phenotype (early onset and bilateral RB). 
However, ICRB staging, an additional clinical feature, did 
not seem to be correlated with genotype. Group E was the 
most frequent staging at diagnosis in our study, which had 
no remarkable differences between patients with null muta- 
tions and in-frame mutations. We believe that tumor staging 
might be more related to access to care rather than the type 
of mutation. 

Interestingly, cases with negative genetic testing results, 
which were mostly non-hereditary (42 unilateral and 3 bilat- 
eral probands), seemed to present with more serious clinical 
staging at diagnosis, with 16 (33.3%) extraocular extension 
tumors. The finding that advanced tumors tend to occur in 
non-hereditary cases might be explained by two hypotheses. 
First, most of the hereditary cases were bilateral RB, and the 
growth of the retinoblastoma tumors in the two eyes was not 
completely synchronous. The first-discovered eye disease 
in bilateral probands likely contributed to early diagnosis of 
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Table 2. Clinical profile of retinoblastoma probands (age and laterality). 



Groups 



probands (no. and %) mean age (months) P value 



laterality (no. and %) 



bilateral 



unilateral 



P value 



null mutation 31 (34.5%) 

in-frame mutation 9 (10.6%) 

negative genetic testing results 45 (52.9%) 



10.7 
13.5 
27.3 



<0.001* 



30 (96.8%) 
6 (66.7%) 
3 (6.67%) 



1 (3.2%) 
3 (33.3%) 
42 (93.3%) 



0.008** 



total 



85 



19.8 



39 (45.9%) 46 (54.1%) 



*ANOVA was applied (F=22.97); there were significant difference between two groups (null versus in- frame, null versus neg, in- frame 
versus neg). **Chi-Square test was applied between null and in- frame (x2=7.0251); 



another eye disease. Second, the growth processes of reti- 
noblastoma tumors might have certain differences between 
patients with and without germ-line RBI mutation, resulting 
in faster tumor growth in probands with negative genetic 
testing. A recent study reported that a previously unrecog- 
nized subtype of retinoblastoma tumor with intact RBI alleles 
had high-level v-myc avian myelocytomatosis viral oncogene 
neuroblastoma derived homolog (MYCN) oncogene amplifi- 
cation, which might initiate retinoblastoma with intact RBI 
alleles [29]. This suggests that the underlying mechanism 
involved in the initiation or development of retinoblastoma 
tumor could be investigated in future studies. 

Distribution of the RBI gene mutation: Our RBI mutational 
distribution results revealed that RBI mutations were scat- 
tered along the genomic sequence (Figure 1), lacking hot 
spots clustered in the important functional domains — A/B 
pocket regions. Richter et al. reported that 13 of 15 missense 
mutations of the RBI gene gathered in the A/B pocket and 
the intervening spacer region, confirming that in-frame 
mutations favored the critical regions [24]. However, only 
five of the nine missense mutations appeared in the A/B 
pocket domains in our study. All five missense mutations 
that appeared in the A/B pocket domains were identified 
in bilateral probands, while only one of the four missense 
mutations that appeared in other regions were identified in 



bilateral probands, suggesting that in-frame mutations in A/B 
pocket domains could be considered to have more complete 
penetrance. Mutations occurring in different regions along 
the RBI gene still need to be investigated in a study with an 
expanded sample size. 

Incomplete detection rate: We reported a detection rate of 
78.6% in our previous paper by using single DNA sequencing. 
Here, a method combining DNA sequencing and MLPA 
increased the detection rate to 92.3% in bilateral RB, which 
is consistent with an earlier study that reported that the RBI 
mutation detection rate was 92.6%> via a combination of full 
sequencing and deletions/duplication analysis of RBI [30]. 
The remarkable enhancing of the effectiveness of detecting 
gross deletions attributable to MLPA was 20% (8/40)), which 
is slightly higher than in previous studies [31]. However, RBI 
defects can also result from hypermethylation of the promoter 
region and different levels of mosaic mutations, as evidenced 
by reports that highly sensitive allele-specific PCR increased 
the sensitivity from 92.6% to 94.8% for detecting mosaic RBI 
mutations [30]. Detecting methylation and mosaic mutations 
of RBI would be interesting in future studies. 

In conclusion, we developed the molecular detection 
project to optimize the clinical management of all patients 
with retinoblastoma. This study reports our results involving 



Table 3. Clinical profile of retinoblastoma probands (staging). 


Groups 


total no. of eyes 






staging (no. and %) 






A 


B 


C 


D 


E 


Extra- 
ocular 


null mutation 


61 


3 (4.9%) 


8 (13.1%) 


16 (26.2%) 


10 (26.4%) 


20 (32.8%) 


4 (6.6%) 


in-frame mutation 


15 


1 (6.7%) 


1 (6.7%) 


0 


4 (26.7%) 


8 (53.3%) 


1 (6.7%) 


negative genetic testing results 


48 


1 (2.1%) 


2 (4.2%) 


3 (6.3%) 


3 (6.3%) 


23 (47.9%) 


16 (33.3%) 


total 


124 


5 (4.0%) 


11 (8.9%) 


19 (15.3%) 


17 (13.8%) 


51 (41.1%) 


21 (16.9%) 



*Fisher's exact was applied between null and in-frame (p=0.261); Additionally, mutations (including null and in- frame) were compared 
with negative testing. Chi-Square test was applied between the intraocular tumor (including staging A, B, C, D, E) and extra-ocular tumor 
(X2=14.9696, p=0.000): 
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Figure 1. Schematic representation of RBI exonic mutation found in this study. Arabic numbers in the parentheses show the occurrence 
times for each mutation. All gray bars represent null mutations, and the bars represent in-frame mutations. The dark gray regions on the 
RBI gene show the pocket domains. 



31 patients who had been identified with null mutations 
and nine patients identified with in-frame mutations. The 
different groups had different clinical outcomes, and the type 
of mutations was related to age of onset and laterality, but not 
related to staging of the tumor, which has implications for 
genetic counseling. In addition, 20% of the total mutations 
were detected with MLPA, which confirms that combining 
sequencing and MLPA is a convenient, reliable, and powerful 
method for evaluating RBI gene mutations in patients with 
retinoblastoma. 

APPENDIX 1. RESULTS OF GENETIC TESTING 
AND CLINICAL FEATURES COVERING 
LATERALITY, AGE AT DIAGNOSIS, STAGING AT 
DIAGNOSIS. 

To access the data, click or select the words "Appendix 1." 
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